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Based on the refined crystal structure comprised of columns of 3x4 planar blocks of NbOg octahe- 
dra and first principles electronic structure methods, we find that orthorhombic (o)-Nbi 2029 intro¬ 
duces a new class of transition metal oxide. The electronic system consists of a large Nb dimer-based 
localized orbital comparable in size to those in organometallic compounds, yet is tightly bound and 
weakly interacting with itinerant electronic bands. These local moments - a rare occurrence for Nb - 
form one-dimensional spin chains that criss-cross perpendicularly oriented conducting “nanowires.” 
The local moment bandwidth is comparable to what is seen in rare earth compounds with extremely 
localized orbitals. The microscopic origin is traced to the local structure of the NbOe octahedra and 
associated orbital-l-spin ordering. The resulting IDs x IDs anisotropic two dimensional Heisenberg- 
Kondo lattice model (s=spin, c=charge) provides a strongly anisotropic spin-fermion lattice system 
for further study. 

PACS numbers: 71.27.-|-a,71.20.Be,72.15.Qm,75.20.Hr 


I. INTRODUCTION 

Transition metal oxides form some of the most intel¬ 
lectually rich electronic phases, viz. the high tempera¬ 
ture superconducting cuprates, the colossal magnetore¬ 
sistance manganites, and ordering of charge, spin, or¬ 
bital and structural changes, especially in 3 d oxides. Nio- 
bates are 4 d oxides that assume a fascinating sequence 
of crystal structures with assorted properties, extending 
from the cubic one-to-one, formally Nb^+, metal NbO to 
the Nb^^-based wide-gap insulator Nb205, and encom¬ 
passing a sequence of mixed valent compounds such as 
Nbi 2029 , Nb22 0 54, Nb 25 062i Nb470ii6, and so on. 

The 12-29 member, Nbi 2029 , displays local moment 
behavior coexisting with a relatively high density of con¬ 
ducting carrierSfi Just how local moments - very rare 
for Nb - emerge and coexist with itinerant electrons in 
this system has remained an enigma for decades. More 
broadly, systems of local moments embedded in conduct¬ 
ing media form a rich platform for unusual phases, with 
phenomena including Kondo systems, heavy fermion 
metals and superconductors, and still unexplained non- 
Fermi liquid behavior. 

Nbi2029 originally attracted attention's^ due to the 
mixed valent character it implied for Nb: 2 Nb'’’’*' -I- 10 
Nb®“'" balancing the oxygen charge, apparently leaving 
two Nb cations with Ad} configuration. It lies just in¬ 
side the line of conducting compounds in the Nb205_2a; 
system, where coexistence of polarons and bipolarons 
had provided the prevailing picture of its conducting 
behavior 4 i^ The tungsten bronzes, which share the same 
underlying WO3 structural motif, have been found to su¬ 
perconduct up to 4 K when heavily dopedi^si 

Nbi2029 is found in two polymorphs, orthorhombic 
(o-) and monoclinic (m—), both based on stacking of the 
same underlying structural feature (described below) but 


being stacked differently. The magnetic susceptibility of 
each polymorph displays Curie-Weiss susceptibility cor¬ 
responding to one spin-half moment per f.uji TO-Nbi2029 
orders around 12 K, in the sense that the susceptibility 
follows a behavior that peaks before dropping, in a man¬ 
ner that can be fit by the Bonner-Fisher form for a one di¬ 
mensional (ID) Heisenberg antiferromagnet 4 o-Nbi 2029 
does not show such a peak down to 2 K>S 

This first (and still one of few) observation of Nb lo¬ 
cal moment behavior, moreover to be coexisting with Nb 
conductivity,—!^ suggested a more specific picture^ of one 
electron becoming localized and magnetic while the 
other is itinerant and Pauli paramagnetic. While this 
compound is a very bad metal, it is nevertheless also 
a standout transparent conductor, displaying resistivity 
^ 3 mf 2 cm with only 10 % variationl^ from room tem¬ 
perature to 2 K, and it possesses a high carrier densityiS 
of one carrier per Nb thus not requiring tuning via ex¬ 
trinsic doping. Fang et al. recently argued from a density 
functional study using the generalized gradient approxi¬ 
mation (GGA) functional that the monoclinic polymorph 
is an itinerant Stoner magnet— i Our studies provide a 
very different picture of o-Nbi2029, and we contrast the 
two pictures below. 

In this paper we demonstrate that orthorhombic o- 
Nbi2 029 presents an unprecedented type of material, 
based on an unusual large but strongly localized and 
magnetic organometallic-like orbital—^ This compound 
consists at the most basic level of an array of linear 
(ID) Heisenberg magnetic chains along the one direc¬ 
tion, crisscrossed by ID conducting nanowires along a 
perpendicular direction, coupled to a two dimensional 
( 2 D) system by Kondo interaction at each site. This new 
IDxlD Heisenberg-Kondo lattice (HKL) system arises 
from the intricate structure of columns of 3 x 4 planar 
units of NbOe octahedra, with the spin-half dimer Mott 



FIG. 1: (Color online) View of the structure in the c — b 
plane (c horizontal) of o-Nbi 2 029. The designations of the 
sites Nbl, Nb2, Nb6 in the 3x4 unit are indicated. The 
dashed lines outline two of the 3x4 units that lie on the same 
level. Other such units (only partially visible) lie at a height 
^ along a. Note that the (green) Nb5 site is the only one not 
on a ‘surface’ of the 3x4 unit. 


insulating system arising from quantum confinement in¬ 
duced extreme localization on Nb dimers with a specific 
4d orbital orientation, while other Nb sites provide elec¬ 
tronic conduction. These results should stimulate further 
study both experimentally, to identify the very low tem¬ 
perature and high magnetic field phases, and theoreti¬ 
cally to probe a new model system for novel magnetic 
and possibly exotic superconducting states, since related 
2D models have been shown to harbor pairing correla¬ 
tions. 


II. STRUCTURE AND APPROACH 

Nbi 2 029 assumes two polymorphs, monoclinic (m- 
^i9ji3-u6 orthorhombic (o-).— Each is based on 
the same 3x4 block structure of corner-linked NbOg oc- 
tahedra (ReOa-type perovskite) that are stacked along 
the a direction with lattice constant a = 3.832 A, thus 
forming rectangular oo x 3 x 4 nanocolumns (in a, b, c or¬ 
der of the Nb sites). The 3x4 blocks contain six crystal¬ 
lographic sites Nbl, ..., Nb6. The two structures differ 
only in how these columns are repeated in the b — c plane, 
a common motif in polymorphic shear structures. 

We focus here on o-Nbi202g^‘^*^ with 6=20.740 A, 
c=28.890 A, whose crystal structure was only settled by 
use of xray diffraction together with convergent beam 
electron diffraction on single crystals by McQueen et alt^ 
The structure at 200 K, shown in Fig. [U has centered 
orthorhombic space group Cmcm (^ 63) with two f.u. 
per primitive cell. The primitive cell consists of two 
such columns, one of which lies ^ lower/higher than the 
other, so these two NbOe columns are connected edge¬ 
wise rather than corner-linked. These shear boundaries, 
we find, severely limit carrier hopping in the b — c plane. 
Due to the distortion of the Nb-0 bond lengths and an¬ 


gles, the 3x4 unit is ferroelectric in symmetry. The cen- 
tering-breflection symmetry operation however results in 
no net electric polarization in the cell. 

Density functional theory (DFT) calculations were 
carried out initially using the spin-polarized Perdew- 
Burke-Ernzerhof GGA functional?^ implemented in 
an all-electron full-potential code fploA^ The on-site 
Coulomb repulsion U was treated within the the GGA-I-U 
approachA2i^ We present results for [/ = 3 eV, the 
value calculated^i for Nb in perovskite SrNbOa, and 
Hund’s energy = 1 eV; we verified that results are 
insensitive to the specific values. A similar value of U 
has been applied in the DFT plus dynamical mean field 
method to study the dynamical excitation spectrum of 
Lii_a;Nb 02 ?^ which as one of the few strongly correlated 
niobates has become of interest because as a strongly lay¬ 
ered transition metal oxide it becomes superconducting, 
thus displaying several similarities to the high tempera¬ 
ture superconducting cuprates. The Brillouin zone was 
sampled with a regular fine mesh up to 657 irreducible k 
points. 

III. THE THEORETICAL PICTURE 
A. Electronic configuration 

We first address the electronic structure before consid¬ 
eration of magnetism. The system was earlier found to 
be too intricate for Hiickel methods to elucidate<^ The 
atom-projected density of states (PDOS) in Fig.[2]reveals 
a sharp and narrow Nb5-dominated peak, with secondary 
contribution from Nb4 and Nb6, that pins the Fermi level 
Ei? and will be found to be unstable to local moment 
formation. Itinerant band states around Ef are domi¬ 
nated by the Nb2 site, with smaller Nb3 participation 
but none from Nb5. The band structure in Fig. [3] reveals 
a Nb52 dimer bonding band lying within the continuum 
of a-axis dispersive conduction bands, with remarkable 
flat-band dispersion of no more than 75 meV along each 
of the three axes. The antibonding partner lies ^^2.8 
eV higher and is bifurcated by additional interactions. 
From the small dispersion of the flat Nb5 band, we ob¬ 
tain the nearest neighbor hopping amplitudes along the 
three axes: tg ~ —10 meV, tg = 17 meV, tg = 0. (The 
dispersion along d is however not simple nearest neighbor 
type, in spite of the small value.) 

The anomalously flat band arises primarily from the 
Nb5 dyz orbital, which significantly has its lobes lying 
in the b — c plane, perpendicular to the column d di¬ 
rection. This orientation provides only a tiny ddd over¬ 
lap (tg above) for hopping along the short d axis. The 
hopping amplitudes given above are orders of magnitude 
smaller than the on-site Coulomb repulsion that will pro¬ 
duce Mott localization of dyz states on the Nb52 dimer. 
The two neighboring Nb5 dyz orbitals form bonding- 
antibonding combinations by strong ddn coupling, split¬ 
ting the molecular orbitals by 2tg^,n. = 2.8 eV. The Nb52 
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FIG. 2: (Color online) Total and orbital-projected density of states for each of the Nb sites in nonmagnetic Nbi 2 029. Note 
that Nb5 dominates the peak at the Fermi level Ef, while Nb2 dominates the occupied states at and below Ef- 
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FIG. 3: (Color online) The full band region of non-magnetic 
Nbi2029, above the gap separating them from the O 2p bands. 
The red horizontal lines indicate a single, remarkably flat, 
dyz bonding and (split apart) antibonding band pair of the 
Nb52 dimer. Symmetry points are X=(-,0,0), Y=(0, f,0), 

z=(0,0,f). 


dimer bonding combination is partially filled and pins 
Ep. 

The Nb2-derived conducting “wires”, a pair on oppos¬ 
ing columnar faces perpendicular to c and two columns 
per primitive cell, lead to four bands that disperse 
through Ef along (fea,0,0). Coupling between the con¬ 
ducting wires lifts the degeneracy, as can be seen from 
the occupied bands along T-X in Figs. [3] and |H Only 
two bands cross Ef along E-Y. There are four flat, ID 
sheets perpendicular to a, plus some small and probably 


unimportant sheets. These flat sheets, with wavevectors 
fcF~±(0.15 — 0.35)^, imply large susceptibilities and im¬ 
pending intra-columnar ID instabilities at several Fermi 
surface nesting calipers. The several different nesting 
wavevectors that connect the various ID FS sheets will 
however tend to frustrate any specific instability {viz. 
charge or spin density waves, or Peierls distortions). 


B. Magnetic character 

The conventional GGA treatment of electronic ex¬ 
change and correlation processes produces a Slater mag¬ 
netic moment arising from the exchange-split flat band. 
However, in flat-band transition metal oxides such as 
this (though rarely as extreme), explicit attention to on¬ 
site repulsion is important, so we focus on the GGA-I-U 
treatmen t that is often successful in treating strong 
interaction effects. Here U is the Hubbard repulsion on 
the Nb ions supplementing the semi-local GGA exchange- 
correlation potentialJi 

From the correlated treatment, one peculiar feature is 
that the magnetic moment arises from a single electron 
on the Nb52 dimer rather than on a single atom. This 
distinction makes it a quarter-filled band system that re¬ 
sults in a dimer Mott insulating (dMI) subsystem;^ a 
phase not yet reported in otherwise conducting 3D sys¬ 
tems. Dimer Mott insulators are rare, although a trimer 
Mott insulator has been reported.— The Nb52 dimer or¬ 
bital density also the spin density, is displayed in 

Fig. [51 and its bonding character (incorporating the in¬ 
tervening O ion) is evident. 
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FIG. 4: (Color online) Fatband plots of the bands of Nbi2029 
near Ef, without and with the strong interaction ([/). The 
magnetic bands within GGA (left panel) highlight the flat 
Nb5 dyz band; the dispersive bands within ±0.6 eV of E_f are 
strongly Nb2 d^j^-derived. The right panel shows the bands 
from a ferromagnetic GGA±U calculation, with lower Hub¬ 
bard band in red. The upper Hubbard band lies above this 
energy window. The splitting of majority (resp. minority) 
conduction bands (yellow, resp. green), ~40 meV near Ef, 
provides a measure of the Kondo coupling of the spins to the 
itinerant electrons. 
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FIG. 5: (Color online) Nb52-centered spin density isosurface, 
looking down on one layer of the 3x4 unit in the & — c plane, of 
the magnetic electron in Nbi 2 029. Since the state is weakly 
mixed in all three direction, the spin density is also the density 
of the corresponding Nb52 dimer centered molecular orbital 
(j>b. A majority of the density lies on the central Nb52 dimer 
dy^ orbitals and the intervening O Py orbital. 


This dMI occupied state however lies within the con¬ 
tinuum of occupied conduction bands. From the nearly 
invisible band mixing with the conduction bands (see 
Fig. |4]), it is found that band mixing (charge coupling) 
with the conduction electrons is quite small. Then, from 
the splitting of majority and minority conduction bands 
(see Fig. 0]) the exchange (Kondo) coupling along b, 
JkSj ■ Sj, between itinerant Sj and localized Sj spins 
at site j, has coupling strength Jx = 20 meV that repro¬ 
duces the spin splitting of the conduction bands in Fig.|3] 
for ferromagnetic alignment. 

Three features may be contributing to the special be¬ 
havior of the Nb5 site. First, quantum confinement: Nb5 
is the only site interior to the 3x4 block, it is not an in- 



FIG. 6: (Color online) Schematic representation of the ac¬ 
tive states in Nbi 2 029. Spin chains (arrows) form columns 
along d but are coupled most strongly along b (hence shown 
antialigned), while they are uncoupled along c. Each spin is 
positioned between two ‘wires’ conducting along the a axis; 
the green shading gives an idea of the conduction electron 
density spread. 


terface atom. Second, its octahedron is much more reg¬ 
ular than for the other sites, having only a single Nb5-0 
distancei^ greater than 2.00 A. Third, the anisotropic dis¬ 
placement parameter {Un describing displacement from 
the ideal position in the a direction) is six times larger 
than for any other Nb site, reflecting dynamic (or possi¬ 
bly spatial) fluctuations in position. 

Nbi 2029 thus presents a novel type of system with in¬ 
teracting local moments in a conducting background - a 
Heisenberg-Kondo lattice system. The earlier ID HKL 
model has been generalized to two-dimensions in a max¬ 
imally anisotropic matter: spin coupling in one direction 
and fermion hopping in the perpendicular direction. The 
basic picture is illustrated schematically in Fig. [6] Geo¬ 
metrically, the picture is one of ID Heisenberg spin chains 
of Nb52 local moments extending along b, sandwiched be¬ 
tween nanowires (or “tapes”) conducting along a, lying 
on the faces of the columns perpendicular to b. The lack 
of dispersion along c can be ascribed both to the elec¬ 
tronic disruption caused by the shear boundaries, and to 
the centering±mirror symmetry operation that serves to 
produce maximum separation of neighboring Nb sites. 

This picture of o-Nbi 2 029 is differs from that of m- 
Nbi 2 029 presented by Fang et aZAi Based on GGA alone 
and the VASP code, they presented a picture of a weak 
itinerant (Stoner) moment arising from a flat band with 
magnetic exchange splitting of 0.3 eV. This result is much 
like what we obtained for o-Nbi2029 at the GGA level 
(see Fig.[4|). Our view is that the flat band requires con¬ 
sideration of correlation effects. Several specific differ¬ 
ences follow, viz. strong local moment leading to a dMI 
subsystem versus a spin-split band Stoner moment in a 
fully itinerant system, and the corresponding coupling to 
the carriers. Such differences can be resolved by more 
extensive experimental study at low temperature. 
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IV. DISCUSSION 

From the hopping amplitudes given above, the largest 
magnetic coupling is not along the structural a chain, 
but rather between neighboring Nb52 pairs in the b 
direction, corresponding to a superexchange coupling 
through the oxide ions at the edges of the 3x4 blocks 
of Jh = 4(t|)^/f/ ^ 0.5 meV = 6 K. We presume the 
direct exchange will be smaller, based on the molecular 
orbital pictured in Fig. 5, because the direct overlap 
with neighboring orbitals seems very minor. The mo¬ 
ments will be relatively uncorrelated along the conduct¬ 
ing a direction (J“ ^ 2 K), with the dominant antifer¬ 
romagnetic (AFM) correlations developing along b only 
at low temperature below J^. This temperature scale is 
roughly consistent with that from susceptibility x(T)^ 
which shows spin-half behavior with 9cw = "14 K, but 
no deviation from Curie-Weiss behavior down to 2 K. 
The monoclinic phase, which we will address elsewhere, 
“orders” via displaying a x peak at 12 K that can be 
fii^ with the Bonner-Fisher form2£ for the ID Heisenberg 
chain (which shows no true long-range order). 

The basic Hamiltonian to model such a system, ne¬ 
glecting band indices and spin chain indices for the sim¬ 
plest model, is 

^ ~ -|- c]j ^Cj-a^a] 

+ Jh -f Jk (1) 

j 3 

where j, a run over the 2D lattice and spin directions, 
a and b are the 2D lattice vectors, a is the Pauli spin 
matrix vector, and a^apCj^p = Sj is the conduc¬ 

tion electron spin at site j in terms of the carrier creation 
operator c]^. The energy scales are t=400 meV, Jh^O-5 
meV, Jif=20 meV. In this model an itinerant electron is 
confined to its nanowire, and suffers only Kondo (spin) 
scattering due to the interaction. The ID Heisenberg 
chains likewise are perturbed only by spin coupling to 
itinerant electrons. 

The strictly ID version of the Heisenberg-Kondo lat¬ 
tice model has been studied actively^liSS in the context 
of high temperature superconducting cuprates, where 
ID structures with entwined charge, spin, and su¬ 
perconducting orders have been detected experimen¬ 
tally. The 2D HKL model has recently been applied to 
model behavior, and pairing symmetry, in heavy fermion 
superconductors.— This new IDcXlDg model of cou¬ 
pled ID charge (c) wires and ID spin (s) chains com¬ 
prises a type of system with frustrating interactions. ID 
conductors have thoroughly studied charge- and spin- 
density wave instabilities at 2kF and Akp respectively. 
Are these instabilities overcome by the Kondo coupling, 
which results in coupled ID conductors and hence a 2D 
system? The AFM ordering tendencies of the Heisen¬ 
berg chain, which never actually attains long range or¬ 
der, is correlated with neighboring chains by the Kondo 


coupling. A plausible mean field picture might be that 
Heisenberg spin chains coupled, either in phase or out 
of phase, with neighboring spins chains by the Kondo 
coupling through the conducting wires. The correspond¬ 
ing spin-exchange and magnetic ordering could intro¬ 
duce superstructure into the fermionic susceptibility at 
commensurate wavevectors, possibly influencing the di¬ 
verging susceptibility at 2kF on the conducting systems, 
through the induced alternating “magnetic field” trans¬ 
ferred through the Kondo coupling. This system may 
lead to a new form of magnetic polarons, accounting for 
the bad conductor behavior that is so prominent in low 
carrier density transition metal oxides. Superconduct¬ 
ing pairing tendencies become an attractive possibility, 
due to the analogies with cuprates; note that Nbi 2029 is 
a magnetic 2D transition metal oxide with '^0.1 carrier 
per metal ion. We expect that low temperature and high 
field studies of this system could reveal unusual trans¬ 
port, thermodynamic, and spectroscopic behavior. 


V. SUMMARY 

In this electronic structure study of o-Nbi2029, two re¬ 
markable features have been uncovered. First, regardless 
of whether strong on-site interaction effects are included 
in the calculation, a remarkably flat Nb 4d band is found 
to exist within a bands of several itinerant (dispersive) 
Nb Ad bands. Secondly, when correlation effects are in¬ 
cluded, the majority spin component of this flat band 
state is fully occupied, but the orbital is shared by two 
Nb sites, the Nb5 sites interior to the 3x4 columnar array 
of NbOe octahedra. The itinerant, conducting states are 
found to reside along the edges of the columns, forming 
ID wires. The magnetic spins are (i) coupled to neigh¬ 
boring spins along a single direction, and (ii) coupled to 
the itinerant carriers by Kondo coupling. This combina¬ 
tion appears to be quite delicate: while the orthorhombic 
polytype that we have studied does not order down to 
2K, the monoclinic polytype with many structural sim¬ 
ilarities “orders” in the sense of a ID Heisenberg chain, 
with the susceptibility following the Bonner-Fisher form. 

The picture of Nbi2029 that arises makes more ex¬ 
plicit the earlier picture of “one magnetic Nb ion and 
one conducting Nb ion.” A local moment does emerge 
as required by the susceptibility data, but it is a highly 
unusual dimer moment rather than a single ion moment. 
Similarly, the conducting state arises not from a single Nb 
ion providing an itinerant electron; instead the conduct¬ 
ing electron is shared between two conducting “wires” on 
either edge of the 3x4 columns of NbOe octahedra. The 
novel state that results is an illustration of how structural 
complexity can give rise to new emergent phases of mat¬ 
ter. This new understanding of this Nb suboxide should 
serve to guide study of the more complex Nb suboxides 
mentioned in the Introduction. 
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